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DEVELOPMENT  OF  A SOIL-WHEEL  INTERACTION  MODEL 
George  Y.  Baladl  (Member,  ISTVS)  and  Behzad  Rohanl 

U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Mississippi 


A3STRACT 


The  development  of  a mathematical  model  for  calculating  the  motion  resis- 
tance, slnkage,  dra%<bar  pull,  torque,  and  side  force  for  a flexible  wheel 
traversing  a yielding  (or  deformable)  surface  Is  described.  In  order  to 
make  the  problem  tractable,  the  deformed  boundary  of  the  wheel  Is  assumed 
to  be  an  arc  of  a larger  circular  wheel.  The  entire  soil-wheel 
Interaction  process  Is  treated  as  two  springs  In  serlos,  one  describing 
the  flexibility  of  the  tire  and  one  doscrlblng  the  elLstlc-plastic 
deformation  of  the  soil.  Mathematical  expressions  are  derived  for  the 
two  spring  constant r in  terms  of  the  load  deflection  characteristics  of 
the  tire,  the  undeflected  configuration  of  the  wheel,  and  the  mechanical 
properties  of  the  soil  (both  shearing  response  and  con5>resslbility 
characteristics) . 

The  systea  of  equations  describing  the  performance  of  the  wheel  is  solved 
numerically  via  a coi^putec  program  called  TIRE.  Using  this  program,  a 
seriev  of  parametric  calculations  is  cotiducted  to  demonstrate  the  applica- 
tion of  the  methodology  and  to  sU^dy  the  performance  of  flexible  wheels 
on  different  types  of  soil  under  various  kinematic  conditions.  A partial 
validation  of  the  proposed  interaction  model  is  established  by  coB^arlng 
the  results  of  a large  nusber  of  laboratory  single  wheel  tests  on  both 
clay  and  sand  with  the  corresponding  model  predictions. 


INTRODUCTION 


The  deterBlnatlon  of  the  response  of  s single  flexible  wheel  traversing  a 
yielding  (or  deforaable)  surface  is  essential  for  the  analysis  of  the 
steering  performance  of  wheeled  vehicles.  Specifically,  the  sinkage, 
motion  resistance,  drawbar  pull,  torque,  and  side  forces  acting  on  a 
powered  flexible  wheel  moving  on  a yielding  soil  must  be  accurately 
determined.  Due  to  the  overwhelming  co^>lexity  of  this  problem,  previous 
research  in  this  area  has  been  directed,  by  and  large,  towards  extensive 
experimentation  and  the  development  of  empirical  equations  relating  the 
various  parameters  of  the  problem^C^ference  1).  Unfortunately,  these 
empirical  equations  ere  not  generic’  and  apply  only  within  the  range  of 
the  experimental  data  on  which  they  are  based.  On  the  other  hand,  most 
of  the  analytical  investigations  conducted  in  this  area  are  based  on  the 
astAimption  of  a rigid  wheel  (Reference  2).  That  is,  the  effect  of  the 
flexibility  (elasticity)  of  the  tire  on  the  kinematics  of  the  wheel  is 
neglected.  Even  in  the  case  of  the  rigid  wheel,  there  is  no  general 
equation  that  can  predict  accurately  the  sinkage  as  a function  of  applied 
load,  configuration  of  the  wheel,  and  the  engineering  properties  of  soil 
(Reference  3).  In  a recent  article,  Fujimoto  (Reference  4)  introduced 
the  flexibility  of  the  tire  in  his  analysis  of  the  performance  of  elastic 
wheels  on  cohesive  soils.  He  introduced  an  emplricsl  relation  between 
the  central  angle  of  the  wheel,  the  internal  pressure  of  the  tire,  and 
tltz  radial  stress  acting  on  the  periphery  of  the  tire.  The  radial  stress 
was  assumed  to  be  constant  over  the  periphery  of  the  tire.  Fujimoto 
concluded  that  the  determination  of  the  radial  stress  is  the  most  diffi- 
cult problem  in  the  analysis  of  soil-wheel  interaction  and  recomsended  an 

empirical  relation  between  the  mobility  cone  index  (Cl)  and  the  radial 
stress. 

The  objective  of  the  present  investigation  is  to  develop  a rational  soil- 
wheel  Interaction  model  that  is  free  from  excessive  emplricisu  and  is 
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general  enough  to  treat  a wide  range  of  problems.  Ihe  core  of  the  model 
Is  'ji  method  for  predicting  the  sinkage  as  a function  of  applied  load, 
deflection  of  the  tire,  slip,  undeformed  geometry  of  the  wheel,  and  the 
lundamental  engineering  properties  of  the  soil  (such  us  cohesion,  angle  of 
internal  friction,  density,  compressibility,  etc.).  Accordingly,  the 
model  can  be  used  to  predict  sinkage  in  sand,  clay,  or  soils  exhibiting 
both  cohesive  and  frictional  properties.  The  equilibrium  conditions  and 
the  sinkage  of  the  wheel  are  then  combined  to  calculate  motion  resis- 
tance, drawbar  pull,  torque,  etc. 


To  demonstrate  the  application  of  the  proposed  model,  a series  of  para- 
metric calculations  is  conducted  to  determine  the  performance  of  flexible 
wheels  on  different  types  of  soil  under  various  kinematic  conditions. 
Also,  a partial  validation  of  the  model  is  established  by  comparing  the 
results  of  a large  number  of  laboratory  single  wheel  tests  on  both  clay 
and  sand  with  corresp'vuding  model  predictions. 

DERIVATION  OF  THE  SOIL-WHEEL  INTERACTIOH  MODEL 
General  Procedure 


The  most  essential  part  of  the  soil-wheel  interaction  model  is  a procedure 
for  determining  the  sinkage  of  a flexible  wheel.  The  basic  parameters 
that  must  be  Included  in  such  a procedure,  are  the  applied  load,  configura- 
tion of  the  wheel,  flexibility  or  elasticity  of  the  tire,  slip,  and  the 
fundamental  engineering  properties  of  the  soil  (such  as  shear  strength 
and  compressibility).  The  development  of  the  physical  soil-wheel 
interaction  model  is  presented  in  detail  in  the  subsequent  sections  and 
is  based  on  the  assumption  that  the  entire  interaction  process  can  be 
simulated  by  two  springs  in  series,  with  one  spring  defining  the  elastici- 
ty of  the  tire  and  the  other  describing  the  elastic-plastic  deformation 
of  the  soil.  These  two  springs  are  then  combined  into  a single  equivalent 
spring  describing  the  interaction  of  the  soil-wheel  system. 

The  sitmilation  of  the  resistance  of  the  soil  by  a spring  constant  leads 
to  a nonuniform  distribution  of  normal  stresses  at  the  soil-wheel 
Interface.  The  shear  stresses  at  the  soil-wheel  interface  are  calculated 
from  a rheological  model  which  describes  the  shearing  stress-strain 
characteristics  of  the  soil.  The  final  step  of  the  analyses  is  to 


determine  the  motion  resistance,  drawbar  pull,  torque,  efficiency,  and 
side  force  for  a flexible  wheel  traversing  a yielding  surface.  These 
parameters  are  calculated  based  on  the  assumption  that  the  deformed 
boundary  of  the  tire  is  an  arc  of  a larger  circular  wheel. 

Spring  Constant  for  a Flexible  Tire 

A typical  load-deflection  curve  for  a flexible  tire  on  a rigid  surface  is 
shown  in  Figure  1 where  A denotes  the  deflection  of  the  tire  at  point  A. 
In  practice,  A is  usually  expressed  as  a percentage  of  the  unloaded 
section  height  of  the  tire  (Figure  2).  The  radial  deflection  of  a generic 
point  B along  the  periphery  of  the  tire  at  an  angle  a is  specified  by 
A^  (Figure  1).  If  Che  deformed  section  of  the  tire  is  characterized  by  a 
continuous  spring  with  constant  , Chen  the  vertical  differential  force 
dF  applied  at  point  B can  be  expressed  as 

dF  • A cosa  da (1) 

t a 


From  Figure  1,  A^  can  be  expressed  in  terms  of  A , o , and  the  unde- 
flected radius  of  Che  wheel  R 

A - R - [cosa  - (l  - (2) 

a cosa  cosa  \ R/J  ' ' 

Substitution  of  Equation  2 into  Equation  1 leads  to 

dF  - Rk^jcosa  ~ (l  - ^)j  da  . (3) 

Also,  from  Figure  1, 

®t  , A 

cos  2-  « 1 - ^ (4) 


In  view  of  Equations  3 and  4 and  static  equilibrium,  the  applied  load  W 
can  be  expressed  as 


Integration  of  Equation  5 leads  to  the  following  relation  for  the  spring 
constant 


‘t  "-7  0^ 

2Rlsin  2 cos  y) 


(6) 


The  spring  constant  can  also  be  expressed  in  terms  of  A by  com- 

bining Equations  4 and  6: 


^ ■ 


2A  Vf  - 1 - (f  - 1)  (1  - I) 


(7) 


Equation  7 is  portrayed  in  the  top  of  Figure  1. 


Spring  Constant  for  Soil 

Let  0^  be  Che  radial  stress  necessary  to  maintain  a slow  expansion  of  a 
spherical  cavity  in  an  elastic-plastic  medium  from  radius  to  R 

(Figure  3a).  The  radial  stress  0^  is  expressed  analytically  in  terms  of 
the  shear  strength  parameters  and  the  volume  change  characteristics  of 
soil  (Reference  5).  The  resistance  of  the  soil  to  expansion  of  the  spheri- 
cal cavity  can  be  simulated  by  a continuous  spring  characterized  by  spring 
constant  k . From  Figure  3a,  the  spring  constant  k can  be  expressed 
as 


where  R - R^  corresponds  to  spring  deflection.  Now  consider  a wheel  of 
radius  R embedded  in  soil  to  a depth  R - (Figure  3b).  The  normal 
stress  at  point  A resisting  the  embedment  of  Che  wheel  is  asstimed  to  be 
equal  to  Che  radial  stress  o^  Inside  the  ej^andlng  cavity.  Similar  to 
expansion  of  the  spherical  cavity  (Figure  3a),  the  resistance  of  the  soil 
to  the  embedment  of  the  wheel  can  also  be  simulated  by  a continuous  spring 
with  constant  k^  given  by 


. ANALOGY  BETWEEN  A WHEEL  EHBEOOEO  IN  SOIL  AND  CAVITY  EXPANSION  PROBLEM 
Figure  3.  Proposed  model  for  computing  the  spring  constant  for  soli 
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k » -r — 0 (9} 

s R - R c 
o 

where  D is  the  unloaded  section  width  of  the  wheel  (Figure  2).  Cond>ining 
Equations  8 and  9 we  obtain 

nfk  + R ) (r  - R ) 

^ — - 1 (10) 

where,  from  Figure  3b 

R + R^  - R^l  + cos  (11) 

R - R^  “ R^l  - cos  (12) 

Substituting  Equations  11  and  12  into  Equation  10  and  solving  for 

cos  6/2  and  6 , we  obtain 

8 8 

co»r‘>/i-?R  

®8  ■ 2 cos'^^l  - ^ (U) 


Substitution  of  Equations  11  and  13  into  Equation  S leads  to  the  following 
expression  for  the  spring  constant  k^  ; 


(15) 


It  is  clear  from  Equation  15  that  the  apparent  spring  constant  of  the  soil 
Is  a function  of  the  engineering  properties  of  soil  through  and  the 

geometry  of  the  tire. 

Equivalent  Spring  Constant  for  the  Soil-Tire  System 

The  model  of  the  soil-tire  system  in  terms  of  the  spring  constants  k^,  and 
k Is  portrayed  in  Figure  4.  The  equivalent  spring  constant  k for 


Based  on  the  concept  of  the  spring  analogy  advanced  In  the  previous  sec- 
tions, the  expression  for  differential  vertical  force  at  a generic  point 
at  the  soil-tire  interface  can  be  expressed  as  (Figure  S) 


Figure  5.  Normal  stress  distribution  along  the 
soil-tire  interface 


Solving  Equation  17  for  o„  , we  obtain 

N 


In  view  of  Equations  9 and  12,  Equation  18  becomes 
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cosa  - cos 



1 - cos  ^ I cosa 

Equation  19  describes  the  distribution  of  normal  stress  at  the  soil- 

tire  interface.  Note  that  at  point  A (Figure  5)  where  a " 0,  Equation  19 
indicates  that  o. o , which  Is  consistent  with  the  assumption  made  in 
the  previous  sections.  On  the  other  hand,  at  the  free  surface  where 
a - tO/2  (Figure  5)  Equation  19  indicates  that  « 0 at  these  points. 


1 

(19) 


Consider  now  a tire  with  turn  angle  n with  respect  to  the  direction  of 
motion.  The  plan  view  of  the  tire  is  shown  in  Figure  6a.  If  slip  in  the 
plane  of  the  wheel  is  defined  by  the  slip  ratio  S , then  slip  in  the 
direction  of  the  motion  can  be  expressed  as 


S 

m 


_S 

co.«n 


(20) 


DIRECTION  OF  MOTION 


MR 

a.  PLAN  VIEW 


« 


■SF 


^N 


b.  STRESS  DISTRIBUTION  ALONG 
SECTION  AA 


Figure  6.  Geometry  of  the  tire  with  turn  angle  n 


The  components  of  shear  stress  parallel  and  perpendicular  to  the  plane  of 


the  wheel  t and 
P 

soil  model  presented  in  Reference  S. 


, respectively,  can  be  obtained  from  the  rheological 
N 


; Vi.*’”..,’  ■ 


Cogibinlng  Equation  20  with  the  rheological  soil  nodel  results  in  the 
following  expressions  for  and 


G| 

[c  + 0^  tan^) 

IS 

GS 

+ C + 

tan^ 

cosn 

G| 

fc  + o„  tani> 

^ M / 

|S  tann 

GS 

+ C + o„  tan^ 
N 

cosn 

where  is  given  by  Equation  19.  In  Equations  21  and  22,  G , C 

II 

4’  correspond,  respectively,  to  shear  modulus,  cohesion,  and  angle  of 
Internal  friction  of  the  material. 


Deflection  and  Slnkage  of  a Flexible  Tire 

If  the  deflection  of  a flexible  tire  on  a rigid  surface  under  a given  load 
W is  denoted  by  A (Figure  1),  then  the  corresponding  deflection  on  a 
yielding  soil  (Figure  7b)  can  be  determined  from  the  concept  of  the 

equivalent  spring  constant 


k + k. 


(23) 


Similarly,  if  is  the  slnkage  of  a rigid  wheel  under  a given  load  W 

(Figure  7c),  then  the  corresponding  slnkage  Z of  a flexible  wheel 
(Figure  7b)  is 


(24) 


The  slnkage  Z^  can  be  calculated  from  the  balance  of  forces  in  the 
vertical  direction  (Figure  8a) 


da 


W - OR 


02  = 0 
2 = 0 


OBIS 


a.  RIGID  SJRFACE-FLEXIBLE  TIRE 


02  < 01 
Z > 0 

Iw 


Z«SINKA6E 


b.  SOFT  GROUND-FLEXIBLE  TIRE 

f 

) 

\ 

SINKAGE  V 

c.  RIGID  TIRE-SOFT  GROUND 


Figure  7.  Variation  of  the  central  angles  0j  and  62 
and  sinkage  Z with  relative  rigidity  of  the 
tire  and  soil 


Hie  shear  stress  : in  Equation  25  can  ba  obtained  from  the  rheological 

soil  model  described  in  Reference  5 and  has  the  following  form 


G^C  + Ojj  taniijs 

icsi  + C + tan^i 


(27) 


where  o„  is  given  by  Equation  26.  The  solution  of  Equation  25  leads  tc 
an  expression  for  6^^  . The  actual  sinkage  can  then  be  calculated 

from  {Figure  7c) 

Z^  « R^l  - cos  . (23) 


Relatlonshlpa  Governing  Single  Wheel  Performance 
Geometry  of  the  problem 

Consider  the  geometry  and  boundary  conditions  for  a flexible  wheel-soil 
system  shown  in  Figure  8b.  The  contact  surface  between  the  wheel  and  the 
soil  is  assumed  to  be  an  arc  of  a circle  with  a radius  equal  to  or  larger 
than  the  undeflected  radius  of  the  wheel  (only  in  the  case  of  the  rigid 
wheel  is  the  radius  equal  to  the  undeflected  radius).  The  center  of  this 
circle  O'  is  located  at  the  intersection  of  the  vertical  line  through 
point  A and  the  bisector  of  the  angle  AOB  . According  to  Figure  7b,  the 
relationship  between  the  angle  0^^  , the  sinkage  Z , and  the  deflection 
of  the  tire  is 


(29) 


Also,  from  the  geometry  of  Figure  7b 


(30) 


From  the  geometry  of  figure  8b 


R’ 


1 - cos(0^  - O^j 


(31) 
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Using  Equations  29  and  30  to  eliminate  Z from  Equation  31,  we  obtain  the 
following  relation  for  R*  In  terms  of  R and  the  central  angles  0^^ 
and  02 


R' 


sin 


01  + 02 


sin 


01-02 


(32) 


1 


Equations  29  through  32  completely  define  the  shape  of  the  contact  surface 
between  the  soil  and  the  tire. 

Tire  Internal  motion  resistance 

The  Internal  motion  resistance  (IHR)  of  the  tire  Is  expressed  In  terms  of 
the  deflection  of  the  tire  on  a rigid  surface.  Data  from  a number  of 
experiments  where  IMR  has  been  measured  are  portrayed  in  Figure  9 (Refer- 
ence 6),  which  shows  that  IMR  Increases  rapidly  with  deflection.  The 
, dashed  curves  In  Figure  9 are  approximate  upper  and  lower  bounds  to  the 
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test  data.  Tlio  solid  curve  In  Figure  9 may  be  viewed  as  the  average 
response  and  Is  fitted  with  the  following  mathematical  expression  for 
calculations  of  internal  motion  resistance: 


'»« ■ [‘(if  + «-‘(i)]  To 


Motion  resistance,  drawbar  pull,  torque,  efficiency  and  side  force 

We  can  now  proceed  to  develop  appropriate  eqttations  for  motion  resistance 
(MR),  drawbar  pull  (DBP),  torque  (T),  and  efficiency  (E).  From  Figures  6 
and  8b 


MR  - R 


-^2) 

'D  f sin^a  + ^jda  + IMR  + 

(V«2) 


MF  cosq 


DBP  - R' 


D J Tp  cos(a  + 


da  - MR 


2 


//  R sine.  \ 

\ Sin-V-^  ) 


(1  - S)  (r  - aJ 


where  and  t are  given  by  Equations  19  and  21,  respectively,  with 


0 replaced  by  MF  - R' 


I / sln20  \ 


sinn  . Sinilarly, 


from  Figures  6 and  8 the  side  force  (SF)  is 
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SF  - R’D  I do  + MR  tann  (38) 

.iiil 

2 


where  Is  given  by  Equation  22.  The  above  system  of  equations  pro- 

vldes  a roiqjlete  solution  to  the  performance  of  a flexible  tire  traversing 
a yielding  soil.  A computer  program  called  TIRE  has  been  developed  which 
numerically  solves  the  above  system  of  equations. 


PARAMETRIC  STUDIES  OF  THE  PERFORMANCE  OF 
A FLEXIBLE  WHEEL  ON  A YIELDING  SOIL 


In  this  part,  the  performance  of  a flexible  wheel  on  both  clay  soil  and 
sand  is  parametrically  investigated  (for  n ■ 0).  In  addition,  the  effects 
of  the  unloaded  section  width,  the  deflection  of  the  tire,  and  the  slip 
ratio  on  the  performance  of  the  wheel  are  also  analyzed.  The  radius  of 
the  flexible  wheel  used  for  the  central  case  is  14.1  in.,  its  w)  Ith  is 
8.28  in.,  and  its  carcass  section  height  is  6.33  in.  All  calculations 
were  conducted  for  an  applied  wheel  load  of  1000  lbs.  The  results  of  the 
parameter  study  are  presented  in  the  following  sections. 

Sinkage 

The  results  of  the  calculations  for  assessing  the  effect  of  soil  type, 
slip  ratio,  and  tire  deflection  on  sinkage  are  presented  in  Figures  10 
through  13.  Figures  10  and  12  indicate  that  for  both  clay  soil  and  sand 
sinkage  Increases  with  increasing  slip  ratio.  The  effect  of  tire  deflec- 
tion on  sinkage  is  portrayed  in  Figures  11  and  13  for  clay  and  sand, 
respectively.  As  indicated  in  these  figures,  the  sinkage  decreases  rapidly 
with  increasing  tire  deflections  from  zero  (rigid  wheel)  up  to  approxi- 
mately 40  percent  deflection.  Beyond  40  percent  deflection,  the  rate  of 
decrease  in  sinkage  is  small. 

Motion  Resistance 

The  effects  of  soli  type,  slip  ratio,  and  tire  deflection  on  motion  resis- 
tance are  shown  In  Figures  14  through  17.  Figures  14  and  16  indicate  that 


51 


Flfur*  10.  UUtionthip  bttwtcn  *lak«Bi 
•nd  slip  r«tlo  (or  clay 


rigura  11.  aaUtlonahlp  batwaan  atnk..-ta 
and  tlra  daflacclon  (or  clay 


rigura  U.  lalatloathlp  batwaan  tlnkaga 
and  altp  ratio  (or  cand 


rifura  13.  RaUtloathlp  batwaan  alokaia 
and  tlra  dadactlon  (or  aand 


Figure  1*.  galatlonahlp  between  notion 
realatance  and  allp  ratio  (or  clay 


Figure  13.  galatlonahlp  batwaan  notion 
realatance  and  tlra  dadactlon  (or  clay 


i 


motion  resistance  initially  decreases  with  Increasing  slip  ratio  up  to  a 
slip  ratio  of  approximately  4 percent  and  increases  thereafter.  This 
initial  decrease  in  motion  resistance  has  been  observed  experimentally  and 
is  attributed  to  the  plowing  action  of  the  tire.  The  Increase  in  motion 
resistance  at  higher  slip  ratios  is  due  to  an  increase  in  sinkage  (see 
Figures  10  and  12).  Relationships  between  motion  resistance  and  tire 
deflection  for  each  soil  type  studied  are  shown  in  Figures  15  and  17.  The 
motion  resistance  initially  decreases  with  increasing  tire  deflection  and 
reaches  a minimum  value  at  about  30  percent  deflection.  At  tire  deflec- 
tions higher  than  30  percent,  the  motion  resistance  increases  again.  The 
Initial  decrease  in  motion  resistance  can  be  attributed  to  the  initial 
rapid  decrei.se  in  sinkage  (see  Figures  11  and  13).  The  increase  in  motion 
resistance  at  deflections  larger  than  40  percent  is  due  to  a rapid  increase 
in  the  internal  motion  resistance  of  the  tire  (see  Figure  9). 

Drawbar  Pull 


Figures  18  through  21  portray  the  effects  of  soil  type,  slip  ratio,  and 
tire  deflection  on  drawbar  pull.  Figure  18  indicates  that  for  clay  soil 
the  drawbar  pull  Increases  rapidly  for  slip  ratios  between  zero  and  about 
10  percent.  For  higher  slip  ratios,  the  increase  in  drawbar  pull  is 
relatively  small.  For  sand,  on  the  other  hand,  the  drawbar  pull  increases 
rapidly  and  reache.s  a peak  value  at  about  20  percent  slip  ratio  (Figure 
20) . The  drawbar  pull  then  drops  for  slip  ratios  in  the  range  of  about  20 
to  50  percent.  Beyond  50  percent  slip  ratio,  the  drawbar  pull  increases 
very  slowly.  This  type  of  behavior  also  has  been  observed  experimentally. 
Relationships  between  drawbar  pull  and  tire  deflection  for  each  type  of 
soil  studied  are  presented  in  Figures  19  and  21.  As  indicated  in  Figures 
19  and  21,  the  drawbar  pull  initially  Increases  with  deflection  up  to  a 
deflection  of  approximately  50  percent.  Beyond  this  deflection,  the 
drawbar  pull  decreases  because  of  a rapid  increase  in  the  internal  motion 
resistance  of  the  tire  (see  Figure  9). 

Effect  of  Section  Width  on  Tire  Performanrg 

Figures  22  through  25  present  the  effect  ot  the  unloaded  section  width  on 
sinkage,  motion  resistance,  drawbar  pull,  and  torque,  respectively,  for 
clay  soil  at  15  percent  tire  deflection.  Figure  22  shows  that  sinkage 
decreases  rapidly  as  tire  width  increases  from  approximately  D/R  “ 0.2  to 
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D/R  “ 0.5.  For  larger  tire  widths  the  decrease  in  ainlcage  is  relatively 
small.  Figure  23  shows  that  the  motion  resistance  decreases  as  the  width 
of  the  tire  increases.  This  is  expected  because  as  the  width  of  the  tire 
increases,  the  sinkage  decreases  (see  Figure  22).  It  should  be  pointed 
out  that  in  Figure  23  the  Internal  motion  resistance  of  the  tire  was 
assumed  to  be  Independent  of  the  width  of  the  tire.  If  the  effect  of 
width  on  the  internal  motion  resistance  of  the  tire  were  taken  into  con- 
sideration, the  result  in  Figure  23  would  have  been  different. 

Figure  24  indicates  that  the  drawbar  pull  increases  as  the  tire  width 
increases.  Most  of  the  increase  in  the  drawbar  pull  takes  place  for  the 
tire  widths  less  than  50  percent  of  the  radius.  For  larger  tire  widths, 
the  rate  of  increase  in  drawbar  pull  is  relatively  small.  This  behavior 
is  also  related  to  sinkage  (Figure  22),  where  it  is  observed  that  most  of 
the  decrease  in  sinkage  takes  place  for  tire  widths  less  than  50  percent 
of  the  radius.  The  relationship  between  torque  and  tire  width  is  shown  in 
Figure  25.  The  trend  in  Figure  25  is  similar  to  Figure  24. 

CORRELATION  OF  TEST  DATA  WITH  MODEL  PREDICTIONS 


Background 

The  results  of  the  extensive  parameter  studies  presented  in  the  previous 
section  Indicated  that  the  model  predictions  are  qualitatively  In  agreement 
with  the  observed  performance  of  flexible  wheels  on  a yielding  soil,  A 
detailed  quantitative  validation  of  the  proposed  model  requires  controlled 
laboratory  tests  and  the  measurement  of  the  appropriate  soli  properties 
discussed  in  Reference  5.  A partial  valldati''n  of  the  model,  however,  can 
be  accomplished  by  using  teat  data  already  documented  in  the  literature. 

The  main  drawback  in  using  existing  data  from  the  literature  is  the  lack 
of  Information  on  the  mechanical  properties  of  the  soil  used  in  the  experi- 
ment. Usually  the  soil  Is  characterized  in  terms  of  simple  indices  such 
as  the  mobility  cone  index  (Cl).  These  indices  must  be  translated  to  the 
appropriate  soli  properties  required  by  the  proposed  model.  This  Is  not 
an  easy  cask  and  requires  a separate  analysis  (divorced  from  the  soil- 
wheel  Interaction  model)  to  make  such  a translation.  Usually  one  is 
forced  to  determine  the  numerical  values  of  several  material  constants 
from  an  index  such  as  the  Cl.  This  inherently  Introduces  uncertainties 
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(or  a bias)  in  the  nunerlcal  values  of  the  consLonts  which,  of  course, 
will  affect  the  degree  of  correlation  between  the  tnodel  predictions  and 
the  test  data.  In  spite  of  such  uncertainties,  a partial  validation  of 
the  proposed  soil-wheel  interaction  model  is  attempted  for  the  zero  turn 
angle. 

Test  Parameters 

Test  data  for  13  different  tires  and  2 soil  types  (clay  and  sand)  were 
selected  from  the  literature  for  correlation  with  model  predictions 
(Reference  1).  The  0/R  of  the  test  tires  ranged  from  0.122  (bicycle  tire) 
to  1.737.  A total  of  163  data  points  was  selected  (65  test  data  for  clay 
and  100  for  sand)  for  different  wheel  loads  and  tire  deflections.  The 
tests,  however,  were  all  conducted  at  20  percent  slip.  Soil  data  for  all 
the  tests  were  given  in  terms  of  the  mobility  cone  index  (Cl).  Using  a 
methodology  developed  in  Reference  7,  the  appropriate  soil  properties 
required  by  the  model  were  estimated  from  the  Cl  data.  A summary  of  all 
the  test  data  and  the  companion  soil  properties  are  given  in  Reference  5 
and  for  the  sake  of  brevity  are  not  included  in  this  paper. 

Model  Predictions 


The  results  of  model  predictions  are  plotted  against  the  corresponding  test 
data  in  Figures  26  through  31  for  sinkage,  drawbar  pull,  and  torque.  Each 
figure  contains  a AS-degree  line  (line  of  perfect  correlation),  a line  of 
least  square  fit,  and  the  standard  deviation  7 which  signifies  the 
deviation  between  the  experimental  data  and  the  corresponding  model  predic- 
tions. It  is  a measure  of  the  deviation  of  the  data  points  in  the  figures 
from  the  line  of  perfect  correlation.  Comparisons  between  the  least  square 
lines  and  the  A3-degree  lines  indicate  that  the  overall  correlation  of  the 
model  predictions  with  the  test  data  is  very  reasonable  in  spite  of  two 
possible  sources  of  error— that  is,  the  general  scatter  in  the  test  data 
and  the  uncertainty  in  estimating  several  soil  properties  from  a single 
cone  index.  The  sinkage,  which  is  one  of  the  most  difficult  parameters  to 
predict,  has  the  lowest  standard  deviation.  The  degree  of  correlation  ! 

exhibited  between  the  test  results  and  model  predictions  indicates  chat  I 

the  physical  basis  of  the  proposed  soil-wheel  interaction  model  is  sound  I 

for  both  cohesive  soils  and  granular  materials.  Therefore,  it  may  be  I 


f concluded  that  the  proposed  model  is  capable  of  simulating  the  Interaction 
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between  a flexible  tire  and  a soil  exhibiting  both  cohesive  and  frictional 
properties. 


SUMMARY  AND  COSCLUSIOMS 

A matheiaatlcal  model  for  calculating  the  motion  resistance,  sinkage,  draw- 
bar pull,  torque,  and  aide  force  of  a flexible  wheel  traversing  a yielding 
soil  has  been  developed  and  computerized  for  numerical  application.  The 
entire  soil-v;heel  interaction  process  was  treated  as  two  springs  In  series, 
one  describing  the  flexibility  of  the  tire  and  the  other  describing  the 
strength  of  the  soli.  Mathematical  expressions  were  derived  for  the  two 
spring  constants  in  terms  of  the  load-deflection  characteristics  of  the 
tire,  the  undeflected  configuration  of  the  wheel,  and  the  mechanical 
properties  of  the  soil.  The  motion  resistance,  drawbar  pull,  torque, 
efficiency,  and  side  force  for  the  flexible  wheel  were  obtained  from  the 
equilibrium  equations  by  assuming  that  the  deformed  boundary  of  the  tire 
is  an  arc  of  a circle  with  a radius  equal  to  or  greater  than  the  unde- 
flected radius  of  the  wheel.  The  model  Is  partially  validated  by  comparing 
the  results  of  a large  number  of  laboratory  test  data  for  single  tires  on 
both  clay  and  sand  with  the  corresponding  model  predictions.  Efforts  are 
presently  underway  at  WES  to  couple  the  soil-wheel  Interaction  model  with 
the  dynamic  equilibrium  equations  of  multi-axle  wheeled  vehicles  for 
analysis  of  the  steering  performance  of  such  vehicles. 
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NOTATION 

Cohesion 

Unloaded  section  width  of  the  tire 

Drawbar  pull  applied  on  the  tire 

Vertical  differential  force 

Efficiency  of  the  tire 

Shear  modulus 

Unloaded  section  height 

Internal  motion  resistance  of  the  tire 

Equivalent  spring  constant  for  soil-tire  system 

Spring  constant  of  the  soil 

Spring  constant  of  the  tire 

Motion  resistance  in  th«*  direction  of  motion 
Motion  resistance  in  the  plane  of  the  tire 
Radius  of  the  tire 


' ^ 

“nr 

R 
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P.^  Initial  radius  of  an  expanded  cavity 

R’  Radius  of  a circle  containing  the  deflected  portion  of  the  wheel 

S Slip  of  the  wheel  in  the  plane  of  the  wheel 

SF  Side  force  Applied  on  the  tire 

Slip  of  the  wha<tl  in  the  direction  of  Motion 

T Torque  applied  on  the  tire 

W Tire  load 

Z Sinkage  of  a flexible  wheel 

Z^  Sinkage  of  a rigid  wheel 

a Generic  angle 

A MaxiMuM  deflection  of  the  tire  on  a hard  surface 

A^  Maxlaua  deflection  of  the  tire  on  a yielding  soil 

A^  Deflection  of  the  tire  at  the  generic  point 

n Angle  between  the  direction  of  notion  and  the  plane  of  the  wheel 


o Radial  stress  inside  a cavity 

c ' 

0.,  Noraal  stress  at  the  soil-tire  interface 

T Shear  stress  at  the  soil-tire  interface 

Shear  stress  perpendicular  to  the  plane  of  the  wheel 

Shear  stress  in  the  plane  of  the  wheel 

4 Angle  of  internal  friction 


